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On the transition from two-dimensional 
to three-dimensional MHD turbulence 

By A. Thess f and Oleg Zikanov J 


We report a theoretical investigation of the robustness of two-dimensional inviscid 
MHD flows at low magnetic Reynolds numbers with respect to three-dimensional pertur- 
bations. We analyze three model problems, namely the flow in the interior of a triaxial 
ellipsoid, an unbounded vortex with elliptical streamlines, and a vortex sheet parallel to 
the magnetic field. We demonstrate that motion perpendicular to the magnetic field with 
elliptical streamlines becomes unstable with respect to the elliptical instability once the 
velocity has reached a critical magnitude whose value tends to zero as the eccentricity of 
the streamlines becomes large. Furthermore, vortex sheets parallel to the magnetic field, 
which are unstable for any velocity and any magnetic field, are found to emit eddies 
with vorticity perpendicular to the magnetic field and with an aspect ratio proportional 
to A 1 / 2 . The results suggest that purely two-dimensional motion without Joule energy 
dissipation is a singular type of flow which does not represent the asymptotic behavior 
of three-dimensional MHD turbulence in the limit of infinitely strong magnetic fields. 


1. Introduction 


Turbulent flows of liquid metals influenced by magnetic fields occur under a wide range 
of circumstances, ranging from metallurgy (Davidson 1999) and fundamental turbulence 
research (Tsinober 2001; Moresco & Alboussiere 2004) to the movement of the Earth’s 
inner core (Moffatt 1978). It is widely believed that when the magnetic Reynolds number 
Rm = RoaUL is small, and the magnetic field is sufficiently strong, homogeneous MHD- 
turbulence becomes purely two-dimensional and the electromagnetic dissipation of kinetic 
energy vanishes. The purpose of the present work is to demonstrate that this belief often 
contradicts to reality, and that homogeneous MHD-turbulence may differ strongly from 
purely two-dimensional turbulence even when the magnetic field tends to infinity. 

When an incompressible fluid with density p, kinematic viscosity v, electrical conduc- 
tivity <7, and permeability po moves with velocity U and characteristic length scale L in 
a uniform magnetic field B , and when Rm <C 1, the state of the flow is characterized by 
the Reynolds number and the electromagnetic interaction parameter defined as 


Re = 


UL 


c jB 2 

~w 


( 1 . 1 ) 


We are interested in fully developed turbulence {Re — >■ oo) far away from walls under 
a strong magnetic field N >> 1. It has been established by Moffatt (1967) (see also 
Davidson 1997) and confirmed experimentally by Alemany et al (1979) and numerically 
by Schumann (1976), Zikanov & Thess (1998), Knaepen & Moin (2004) that the turbulent 
flow tends to become two-dimensional so as to avoid electromagnetic (Joule) dissipation. 
Sommeria & Moreau (1982), who identified the importance of walls perpendicular to the 
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Figure 1. Wavenumber distribution of kinetic energy at N >> 1. 


magnetic field, have aptly summarized our conceptual understanding of MHD turbulence 
at N > > 1 by characterizing the distribution of kinetic energy among the wavenumbers 
k± and perpendicular and parallel to the magnetic field, respectively. This distribution 
of kinetic energy in the state of so called “quasi-two-dimensional turbulence” is sketched 
in figure 1. For N >> 1 the kinetic energy of the flow is confined to a narrow region in the 
wavenumber space, outside the “Joule cone” in which electromagnetic energy dissipation 
takes place and below those wavenumbers k± for which viscous dissipation dominates. 
The angle (p ~ iW 1 / 2 of the energy containing region shrinks to zero as N — > oo. This 
picture has tempted researchers to assume that for sufficiently strong magnetic field 
quasi-two-dimensional MHD turbulence becomes purely two-dimensional and the Joule 
energy dissipation vanishes. 

However, there have been strong experimental and numerical indications which sug- 
gest that this hypothesis is unlikely to be correct. In their experiments on freely decaying 
homogeneous MHD turbulence in mercury Alemany et al. (1979) and Caperan & Ale- 
many (1985) found that velocity fluctuations parallel to the applied magnetic field and 
strong Joule dissipation persisted even at high interaction parameters. Disagreements 
between drag measurements of MHD flows in channels with very large aspect ratios (low 
influence of Hartmann walls) (Tsinober 1990) and purely two-dimensional numerical 
simulations of turbulent channel flow by Jimenez (1990) further support the view that 
quasi-two-dimensional turbulence may never become purely-two-dimensional. Direct nu- 
merical simulations of forced MHD turbulence by Zikanov & Thess (1998) indicated that 
two-dimensional vortices, once they have formed, may undergo violent three-dimensional 
instabilities resulting in an intermittent behavior. Finally, numerical simulations of MHD 
turbulence subject to a two-dimensional forcing, which have been carried out by Nakauchi 
et al. (1992) using an EDQNM model, have provided evidence that even for very high 
values of the interaction parameter and strong non-isotropy of the flow, the Joule energy 
dissipation of the flow tends to a nonzero finite value. 

The foregoing observations lead us to the conclusion that purely two-dimensional tur- 
bulence with zero Joule dissipation is a singular state which may never be reached in 
MHD, similarly as ordinary hydrodynamic turbulence in the limit of vanishing viscosity 
will never come close to a state of zero viscous dissipation. We may formulate the follow- 
ing hypothesis: When an electrically conducting fluid is in the state of fully developed 
turbulence (u — » 0) the Joule dissipation of kinetic energy will remain finite as B — >■ oo. 

In this paper we investigate the singular character of purely two-dimensional turbulence 
using stability analysis of model problems. We shall demonstrate that for arbitrarily 
strong magnetic fields a two-dimensional inviscid flow has always access to a sufficient 
number of unstable degrees of freedom to escape from purely two-dimensional behavior. 
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2. A didactical model for the instability of purely two-dimensional MHD flow 

2.1. An exact solution for inviscid MHD flow in a triaxial ellipsoid 

Consider an inviscid electrically conducting incompressible fluid which is confined to the 
interior of the triaxial ellipsoid 


a_ , V_ £2 

a 2 b 2 c 2 


(2.1) 


and subjected to the influence of a homogeneous magnetic field B = Be z . If the magnetic 
Reynolds number is small, the dynamics of the flow is governed by the equations 


tt + (v ■ V)v = --Vp+ -J x B, V • v = 0 (2.2) 

at p p 


J = a(—V(j) + v x B), V • J = 0 (2.3) 

where J is the electric current density and <f> the electric potential (Roberts 1967; Moreau 
1990; Davidson 2001). The equations are supplemented by the boundary conditions at 
impermeable electrically insulating walls of the ellipsoid 

v • n = J • n = 0. (2.4) 

If there is no magnetic field, there exists a family of three-dimensional time-dependent 
exact solutions of the Euler equations 

. [ c b 1 TT , r a ci TTT . \ b a 

v{x,y,z,t) =U{t) y-e z -z-e y + V(t) \z-e x — x-e z \ + W(t) x-e y -y-e x , 

0 C L C Cl -I (X 0 

(2.5) 

where ( U,V,W ) satisfy the nonlinear ordinary differential equations, which are mathe- 
matically identical to the equations describing the free rotation of a solid body (Kerswell 
2002). The velocity field described by (2.5) automatically satisfies the incompressibility 
constraint, the free-slip condition, and is a superposition of three basic two-dimensional 
flows. Each of them has elliptical streamlines and a spatially uniform vorticity directed 
along the axis of the ellipsoid which is perpendicular to the plane of motion. 

We extend the existing theory to the case when a magnetic field is present by expressing 
the electric current density in the form 

,, . T , \ \ c b 1 . . r a ci . . I" b a 1 . . 

J {x,y,z,t) = I{t) y-e z -z-e v + J(t) \z- e x - x- e z \ + K(t) x-e y -y-e x (2.6) 
b c J Lc a J [ a o 

which satisfies the condition V • J = 0 as well as the boundary condition. The coefficients 
I, J, K describe the strength of eddy currents induced by the movement of the fluid in 
the magnetic field. Taking the curl of Ohm’s law (2.3) we obtain, after some algebra, 
I = aBVab/(b 2 + c 2 ), J = -aBUab/(a 2 + c 2 ), K = 0. The vorticity f2 = V x v 
corresponding to the velocity field (2.5) 

n = u{t) {l + l) ex+ m (I + D ey + w{t) (? + ez (2 - 7) 

satisfies exactly the vorticity equation 


^ = (fi-V)v+-Vx(JxB) 

at p 


(2.8) 
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provided that the coefficients U, V, W obey the equations 

(b 2 + c 2 )U = ( b 2 - c 2 )VW - ) U (2.9) 

(c 2 + a 2 )V = ( c 2 - a 2 )UW - ^ V (2.10) 

(a 2 + b 2 )W = (a 2 - b 2 )UV (2.11) 


It should be emphasized that the derivation of this system from the governing equations 
(2.2) and (2.3) did not involve any approximation. 

The system (2.9)-(2.11) can be made nondimensional by introducing dimensionless 
variables according to (U, V, W) = a(C/*, V*, IF* ) and t = t*/a where 
a = (aB 2 / p)a 2 b 2 /[(a 2 + c 2 )(b 2 + c 2 )] is the inverse of the Joule decay time. Using the 
abbreviations A—a/c and B = b/c and dropping the asterisk the equations become 


U = 
V = 

w = 


B 2 


1 


B 2 + 1 
1 -A 2 
1 + A 2 
A 2 -B 2 
A 2 + B 2 


vw-u 


uw-v 


uv 


(2.12) 

(2.13) 

(2.14) 


This nonlinear system has a number of remarkable properties. The nonlinear terms con- 
serve the total kinetic energy and the total angular momentum 


E = i [(1 + B 2 )U 2 + (1 + ^4 2 )U 2 + (,4 2 + B 2 )W 2 } (2.15) 

L=( 1 + B 2 ) 2 U 2 + (1 + *4 2 ) 2 V 2 + (*4 2 + B 2 ) 2 W 2 (2.16) 


The magnetic field gives rise to an anisotropic damping, embodied in the linear dissipative 
terms on the right-hand-side of equations (2.12) and (2.13). The total kinetic energy obeys 


dE 

dt 




(2.17) 


where e = (1 + B 2 )U 2 + (1 + *4 2 )V 2 is the rate of Joule dissipation. As a result, the 
kinetic energy of the flow will always decay unless it is in a purely two-dimensional 
motion perpendicular to the magnetic field (U = V = 0, IF ^ 0). Observe that the 
electromagnetic interaction parameter N (cf. equation (1.1)) does not appear in the 
system because we have used the Joule timescale a -1 to nondimensionalize the velocity. 
Instead, the values of U, V, and W represent the order of magnitude of N. 


2.2. Linear stability of steady states 

Equations (2.12)-(2.14) admit a steady solution (U. V. IF) = (0,0, Wo) describing two- 
dimensional motion perpendicular to the magnetic field. We analyze the linear stability 
of this solution by solving the equations for infinitesimal perturbations (£, f?,C) 


(B 2 - 1 \ 

«=-« + U*n Wo ) 

/ I _ A ' 2 \ 

C = o. 


(2.18) 


V 


(2.19) 

(2.20) 
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The solvability condition yields the growth rates A in the form 

Aj -° (2 - 2i> 

Instability occurs when A>l,.B<lorA<l,.B> 1, i.e. when the magnetic field 
is parallel to the middle axis of the ellipsoid and the IT-component exceeds the critical 
value 

w I (A 2 + 1)(1 + B 2 ) . . dim oB 2 

W c = j — -p- or, m dimensional form, W? = 

y (A 2 - 1) (1 - B z ) p 

(2.22) 

When a — > oo and b -» 0 we have W c — > 1 and IT* m — > aB 2 / p which implies that a very 
elongated sheet-like structure becomes unstable when the interaction parameter defined 
as N = aB 2 /pW is smaller than N c = 1. 

The instability that we have identified in this particular MHD problem is the elliptical 
instability, widely known in ordinary hydrodynamics as a fundamental mechanism for 
the transition to turbulence (see, e.g., Kerswell 2002). It is generated by parametric 
resonance between Kelvin waves, coming from the homogeneously rotating part of the 
velocity field, and the strain field, being a result of the elliptical shape of the streamlines. 
We have just shown that a magnetic field can damp but not completely suppress this 
instability. 

2.3. Nonlinear dynamics and the singular character of purely two-dimensional motion 

We now analyze the nonlinear evolution of the system focusing the attention on the case 
where the initial state is an almost two-dimensional flow characterized by U( 0) <C IT (0) 
and T(0) <C W(0). Solution of (2.12)-(2.14) for A = 2 and B = 0.5 shows that as 
long as IT(0) < W c = 5/3, the system always relaxes to a purely two-dimensional state 
( U = 0, V = 0, IT < IT C ). Figure 2 shows that the evolution from an initial condition 
with IT (0) > IT C proceeds in a quite different way. Although the system finally settles 
at a purely two-dimensional state, the evolution towards this state is characterized by 
long periods of nearly two-dimensional motion, occasionally interrupted by violent three- 
dimensional transients. During these transients kinetic energy is fed from the nondis- 
sipative two-dimensional mode IT to the three-dimensional modes U and V as can be 
seen in figure 2c. The transients lead to significant Joule dissipation of kinetic energy 
as illustrated in figures 2b and 2d. Decay of initially almost two-dimensional MHD flow 
is accompanied by strong divergence from two-dimensionality due to the action of the 
elliptical instability. 

The singular character of purely two-dimensional evolution is further illustrated by 
extending our model so as to include a phenomenological forcing term, which corresponds 
to application of an external torque in the ^-direction. Equation (2.14) is replaced by 

4 2 _ n2 

w = vv + 1 (2 ' 23) 

The new model admits an exact two-dimensional solution 

[7 = 0, V = 0, IT = t (2.24) 

whose kinetic energy grows monotonically with time and whose Joule dissipation is zero. 
However, the numerical solution computed from a weakly three-dimensional initial con- 
dition shows a completely different behavior. The flow corresponding to the numerical 
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Figure 2. Decaying solution of (2.12-2.14). (a), Two-dimensional component W(t.). (b). Total 
kinetic energy (2.15). (c), Three-dimensional components U(t) (solid line) and V(t ) (dashed 
line), (d). Rate of Joule dissipation e. 




Figure 3. As in Figure 2 but for forced system (2.12, 2.13, 2.23). 


solution shown in figure 3 is periodic in time and involves long lasting periods of two- 
dimensional motion interrupted by short three-dimensional excursions. Most importantly, 
the flow is ” statistically” steady in that the time-averaged kinetic energy does not change 
and the mean Joule dissipation is finite. The flow is therefore neither mathematically nor 
physically close to the two-dimensional behavior suggested by the exact solution (2.24). 
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Although our simple model is only a caricature of the behaviour of real MHD turbu- 
lence, it allows us to draw a number of conclusions which serve as a guideline for the 
investigation of the full problem. In particular, we have learned 

• that a purely two-dimensional MHD flow is prone to the elliptical instability, 

• that the nonlinear three-dimensional evolution of the elliptical instability effectively 
extracts kinetic energy from the system on account of Joule dissipation, 

• that an initially nearly two-dimensional flow does not necessarily stay close to purely 
two-dimensional evolution. 


3. Instability in two-dimensional turbulence in a strong magnetic field 

For further analysis it is convenient to rewrite (2.2) and (2.3) as (Roberts 1967) 


d\ 1 oB 2 ,J 2 t 

+ v-V v = — Vp A l — , 

at p p az z 


V ■ v = 0. 


(3.1) 


We are interested in the stability of two-dimensional flows for which the Joule dissipation 
vanishes. The most general flow of this type is of the form 


v(x, V, t) = U(x, y,t)+W (, x , y, t)e z 


(3.2) 


where U has the properties U • e 2 = 0, V • U = 0 and U and W are solutions of 

f + (U-V)U = -ivP, ^ + (U-V)W = 0 (3.3) 

The two terms in (3.2) describe the superposition of a two-dimensional motion perpen- 
dicular to the magnetic field with a motion parallel to the magnetic field. 

Our hypothesis formulated in section 1 implies that the solution (3.2) of equation 
(3.1) is always unstable in the sense that for any initial condition containing a weak 
three-dimensional perturbation the solution will diverge from the two-dimensional one 
and evolve towards a solution with nonzero mean Joule dissipation no matter how strong 
the magnetic field is. A rigorous proof of this assertion is a formidable task, as it would 
require consideration of the initial value problem (3.1) for all admissible fields U and V. 
We shall simplify the task considerably by considering two families of structures which 
we believe are representative of the solutions to (3.3) namely, columnar vortices with 
elliptical streamlines and vortex sheets whose velocity is parallel to the direction of the 
magnetic field. 


3.1. Instability of motion perpendicular to the magnetic field 
We consider an unbounded strained vortex described by 

U (x,y) = — f lEye x + f lE~ l xe y (3.4) 

When the eccentricity E = 1, the flow is a solid body rotation around the axis of the 
magnetic field with vorticity 2D. When E > 1, the flow has a uniform strain and its 
streamlines are ellipses with semi- axes a = \[E and b = 1/a/E- We extend the analysis 
of Bayly (1986) and Waleffe (1990) to the magnetic case by perturbing the basic flow 
(3.4) according to v = U + u. Linearizing (3.1) with respect to the perturbation we obtain 

_ 1 d 2 u 
d z** 


^ + (U ■ V)u + (u ■ V)U = --Vp - aA 

(JX p 


V • u = 0. 


(3.5) 
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where a = aB 2 / p. This system admits solutions of the form 


{u, p} = {u, p}(t) - exp 



(3.6) 


Inserting (3.6) into (3.5) we obtain equations identical to those for the nonmagnetic 
problem solved by Bayly (1986) and Waleffe (1990). 

Uj + ikixiiij + ik m A m iXiUj + Ajiui = -ikjp, kjUj = 0, (3.7) 

where A 12 = —HE, A 21 = fi /E and Aij = 0 otherwise. The solution includes the 
wavenumber vector k rotating around the axis of the vortex as 


k (t) = k 0 [sin 6 cos <j>(t ) , E sin 0 sin <j>(t ) , cos 0] 


(3.8) 


where <f>(t ) = £l(t — to), and the Floquet solution u (t) = e At w [<j>(t)\, where w is a 27r- 
periodic function and the growth rate is 

\(Sl,E,6) = ^ln[p(E,0)}. (3.9) 

Here, p are eigenvalues of an auxiliary eigenvalue problem. The amplitude of perturba- 
tions is 


with 


F(t) = A {E,9,Q)t-a [ 
Jt , 


|u(t)| = |w(t)| ■exp[F(t)] 
cos 2 Odt' 


to sin^ 6 cos 2 Sl(t* — t' 0 ) + E sin^ 6 cos 2 f l(t' — t ! 0 ) + cos 2 6 


(3.10) 

(3.11) 


The second term in the right-hand-side is the only correction to the non-magnetic so- 
lution. Reverting to dimensionless time <f> = fi(i — to) and using a magnetic interaction 
parameter defined by 

N = a/n = aB 2 / P n (3.12) 


we can express F as 


m = 


2tt 


In \p{E,0)\-N 



d£ 

1 + tan 2 6*[cos 2 <f>' + E sin 2 cp 1 ] 


(3.13) 


The system is unstable if F increases over one period, i.e. F(2n) — F( 0) > 0, it is stable 
if F(2n) - F( 0) < 0 and the condition for neutral stability is therefore F(2n) - F( 0) = 0. 
From the last condition we obtain the desired expression for the neutral surface as 


N(E,e) = In \p(E,6)\ • { [ 2 \ , } (3.14) 

{Jo 1 + tan 0[cos z <p + E sm q>\ J 

Figure 4a shows the numerically computed values of (3.14). For E = 1 we have N = 0 
indicating that a vortex with circular stream lines is always stable. The modes take the 
form of Kelvin waves which rotate about the magnetic field with a frequency 2f2 cos 6 (see 
e.g. Greenspan 1968). At E > 1, the flow becomes unstable with respect to perturba- 
tions located within a band 9~{E) <0< 9 + {E). Notice that the location of this unstable 
band is the same as in the nonmagnetic problem since it is determined by the zeros of 
In [p(E, ff)\. In particular, this band originates at 6 — n/3 (Bayly 1986). For a given value 
of E > 1 the maximum of N{E,6) over all 6 in the unstable band defines the critical 
magnetic interaction parameter N C (E) and a critical angle 0 C (E). These quantities are 



2D-3D transition in MHD turbulence 


71 



o L 


1 1.5 2 

2.5 

E 

3 3.5 

(b) 


(c) 

- / 

0.38 


Stable 

0.36 


/ 

/ 

<x> 

0.34 


Unstable 

0.32 



3 4 1 2 

E E 


Figure 4. Stability of an unbounded strained vortex, (a), Neutral magnetic interaction param- 
eter (3.14) as a function of eccentricity E and angle 6 between the wavenumber vector k of the 
perturbations and the axis of the vortex, (b) and (c), critical interaction parameter N c and 6 C 
as functions of E. 


plotted in figure 4b and c. Figure 4b shows that N c increases rapidly with increasing ec- 
centricity while the critical angle 6 C stays in the vicinity of 7t/3. The monotonic increase 
of the critical interaction parameter implies that strongly elliptical (almost sheet-like) 
structures are particularly prone to instability. Having in mind that the nonlinear evo- 
lution of initially circular vortices governed by the two-dimensional Euler equation (3.2) 
always leads to the generation of strongly elongated vortex sheets, we can conclude that 
the two-dimensional evolution proceeds in a way as to make the system more susceptible 
to the elliptical instability. The fact that 0 C ~ 7t/3 may seem counterintuitive, as one 
usually expects in MHD that unstable structures have a tendency to align with the mag- 
netic field ( 6 — » 7t/2). However, it should be noticed that 6 C is monotonically increasing 
and probably asymptotes towards 7t/2. 

Our stability result can be also viewed from another perspective. In a given magnetic 
field any forced two-dimensional flow will eventually become unstable, once the vorticity 
amplitude has reached a sufficiently high level fi such that N falls below its critical value 
and the system becomes three-dimensionally unstable. 

3.2. Instability of motion parallel to the magnetic field 
In purely two-dimensional flow the velocity component parallel to the magnetic field 
behaves like a passive scalar, transported by the velocity U (x,y,t), cf. equation (3.2). 
Such evolution is known to produce steep gradients of W ( x , y, t) [see e.g. Kraichnan & 
Montgomery (1980), Lesieur (1990)]. Since the archetype of such structures is a single 
vortex sheet, which is known to undergo a Kelvin- Helmholtz instability, we will discuss 
the stability of the flow 

W(x) = U 0 


(3.15) 





72 


A. Thess & 0. Zikanov 


under the influence of the magnetic field. This problem, along with more general unidi- 
rectional velocity fields, was originally considered by Drazin (1960). However, the spatial 
structure of the unstable modes were not investigated there and is therefore the central 
focus of the discussion given below. Infinitesimal perturbations u superimposed upon the 
basic flow (3.15) are governed by the linearized equation 


<9u 9u dW 

~7u + f u x — e 2 

at oz dx 


--Vp- a A - 1 
P 


d 2 u 

dz*’ 


V • u = 0. 


(3.16) 


Squire’s theorem, which holds for flow in a parallel magnetic field (Michael 1953; Stuart 
1954), permits one to reduce the analysis to two-dimensional perturbations of the form 
u = (u, 0, w). They can be decomposed into normal modes according to 


u(x, z, t) = u(x) ■ exp [ik(z — ct)\, . . . 


(3.17) 


After introducing a stream function ip via u = ikip and w = —dip/dx, nondimensional 
variables via U = UoU *, x = x*/k, c = Doc*, t = t' * /Uok, an interaction parameter 
N = aB 2 /pUok, and dropping the asterisk, it can be shown (Drazin 1960) that the 
complex wave velocity c is determined by the third-order polynomial equation 


-ic( 1 + c 2 ) + — (1 + 3c 2 ) = 0. 
The normal mode for x > 0 is given by 

ip(x) = exp 



(3.18) 


(3.19) 


where it is assumed that the square root with the positive real part has been taken. 

Drazin (1960) has demonstrated that all three roots of equation (3.18) are purely 
imaginary and that one of the solutions has always positive imaginary part, but he 
did not explicitly calculate c(N) and ip(x). Figure 5a shows the imaginary parts of the 
three roots. The wavenumber k does not explicitly appear in the solution because it 
represents the only characteristic length scale of the problem and has therefore been 
used to nondimensionalize the equations. However, the wavenumber enters the definition 
of N. 

For each value of N there are two stable and one unstable solutions. When N=0 we re- 
cover the classical Kelvin- Helmholtz instability with c — i. For N > 0 the imaginary part 
of the solution originating from c = i always remains positive which implies that the vor- 
tex sheet is unstable for arbitrarily strong magnetic fields. The modes with small N, which 
correspond to small perturbation wavelengths and small magnetic fields are growing 
fastest, while the modes with large N, corresponding to large scales and strong magnetic 
fields grow slowest. The asymptotic behavior of the wave velocity for the unstable mode is 
c = i(l - N/A) for N « 1 and c = i/y/{ 3) for N » 1. The spatial structure of the un- 
stable modes is best revealed by the stream function tp(x. z, t) = Re{ip(x) ■ exp[i(^ — ct)]} 
(nondimensional variables) which, for t = 0, is (see figures 5b and 5c) 


ip(x , z) = exp 





' N 

K 1+ 4, 

) X 

COS 

-—x + z 

4 


ip(x, z) = exp 


_ 3 1 / 4 


(?) 


1/2 


COS 


— 3 3 / 4 


(?) 


1/2 


X + z 


for N <1 (3.20) 
for N > 1 (3.21) 
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(a) 



(b) N=0.1 



(c) N=1 0 



Figure 5. Instability of a vortex sheet, (a), Solid lines - roots of the dispersion equation (3.18). 
Dashed lines - asymptotic limits il/3 1 ^ 2 , —3N/4. (b), Streamfunction of the unstable mode at 
N <C 1 (3.20). (c), Streamfunction of the unstable mode at N ;§> 1 (3.21). 


The unstable modes for large N have an aspect ratio z/x ~ iV 1 / 2 as is characteristic of 
MHD flows in strong magnetic fields. 
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